Abstract. We present a scenario to explain the lithium- 
rich phase which occurs on the red giant branch at the 
so-called bump in the luminosity function. The high trans- 
port coefficients required to enhance the surface lithium 
abundance are obtained in the framework of rotation- 
induced mixing thanks to the impulse of the important 
nuclear energy released in a lithium burning shell. Under 
certain conditions a lithium flash is triggered off. The en- 
hanced mass loss rate due to the temporary increase of 
the stellar luminosity naturally accounts for a dust shell 
_ , formation. 
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1. The lithium-rich RGB stars 

Since its discovery by Wallerstein & Sneden (1982), vari- 
ous scenarii have been proposed to explain the lithium-rich 
giant phenomenon. Some call for external causes, like the 
engulfing of nova ejecta or of a planet (Brown et al. 1989; 
Gratton & D'Antona 1989; Siess & Livio 1999). Others re- 
fer to fresh lithium production by the Cameron & Fowler 
(1971, hereafter CF71) mechanism associated to a trans- 
port process in the deep radiative layers of the red giants 
(Fekel & Balachandran 1993; de la Reza et al. 1997; Sack- 
mann & Boothroyd 1999; Jasniewicz et al. 1999; Charbon- 
nel & Balachandran 2000, hereafter CBOO; Denissenkov & 
Weiss 2000, hereafter DWOO). The latter explanation is 
sustained by the depletion of Be in the atmosphere of the 
lithium-rich giants (Castilho et al. 1999), and by the loca- 
tion of these objects on the Hertzsprung-Russel diagram 
(CBOO). Indeed Hipparcos parallaxes allowed the locali- 
sation of the short episode of lithium production at the 
bump of the luminosity function on the red giant branch 
(RGB). CBOO thus concluded that the lithium- rich RGB 
phase was a precursor to the mixing process which leads 
to other abundance anomalies, like the low carbon isotopic 
ratios seen in ^ 96% of the low-mass stars (<2.2 - 2.5M0; 
Charbonnel & Do Nascimento 1998). This brought very 
tight constraints on the underlying mixing process and 
in particular on its extreme rapidity. It put in difficulty 
models where high lithium abundances can be reached at 
any time between the RGB bump and tip (Sackmann & 
Boothroyd 1999, DWOO). 

In this paper we propose a scenario which accounts 
for the high transport coefficients required to enhance the 
surface lithium abundance of RGB stars in the framework 
of rotation-induced mixing (§2). We show that under cer- 
tain conditions a thermal instability can occur, triggered 
by the nuclear burning of the freshly synthetized lithium 
in a lithium burning shell (LiBS) external to the classical 



Send offprint requests to: Corinne Charbonnel 
corinne.charbonnel@obs-mip.fr 



hydrogen burning shell (HBS; §3) |]. The enhanced mass 
loss related to the increase of the stellar luminosity during 
the lithium-rich phase can explain the circumstellar dust 
shells observed around some of these objects. 

2. The scenario 

2.1. Rotation induced mixing 

Observational evidence has accumulated on the existence 
of an extra-mixing process which becomes efficient on 
the RGB when the low-mass stars reach the luminosity 
function bump (Brown 1987; Gilroy & Brown 1991; Pila- 
chowski et al. 1993; Charbonnel 1994; Charbonnel et al. 
1998; Gratton et al. 2000). At this point the outwardly- 
moving HBS crosses the mean molecular weight disconti- 
nuity created by the first dredge-up. As shown in Charbon- 
nel (1995) and Charbonnel et al. (1998), an extra-mixing 
process related to rotation and previously inhibited by a 
strong molecular weight gradient can then easily connect 
the convective envelope to the deeper layers where the 
nucleosynthesis occurs and consequently lower the carbon 
isotopic ratio. Rotation was also invoked to account for 
C,N,0 and Na anomalies and '^He and ^Li destruction 
in low-mass and low-metallicity RGB stars (Sweigart & 
Mengel 1979; Charbonnel 1995; Denissenkov & Tout 2000, 
hereafter DTOO). 

Recently, DWOO investigated the effect of rotation- 
induced mixing on the surface lithium abundance using 
prescriptions based on Zahn and Maeder developments 
(Zahn 1992, Maeder & Zahn 1998). They used the cor- 
responding transport coefficients Dr derived by DTOO for 
globular cluster low-mass giants. They found that these 

^ Our scenario is different from the previous works which ex- 
plored the possible occurrence of shell flashes on the RGB and 
which considered only the reactions of the CNO cycle (Bolton 
& Eggleton 1973, Dearborn et al. 1975, Von Rudloff et al. 
1988). It also differs from the model developed by Fujimoto 
et al. (1999) to explain Mg and Al anomalies in globular clus- 
ter red giants. These authors indeed proposed that a flash is 
triggered off deeper in the HBS due to the inward mixing of H 
down into the degenerate He core. 
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values of Dr were too low to lead to lithium enrichment 
and invoked the engulfing of a planet in order to trigger 
and dope the transport process. While this combined sce- 
nario is rather attractive, it is harmed by the fact that it 
can happen at any time on the RGB while lithium-rich gi- 
ants with solar metallicity discovered up to now do stand 
by the bump (CBOO). However what was not investigated 
by DTOO and DWOO is the structural response of the star 
to the extra-mixing. Their computations of the transport 
coefficients and of the abundance variations are indeed 
done in a post-processing approach where standard stel- 
lar structures are used as background models. We propose 
here that the structural response to the mixing is actually 
the cause of the increase of which is necessary to en- 
hance the surface lithium abundance. 

2.2. Impact of the nuclear energy released by lithium 
burning 

Fig. |l] shows the abundance profiles of the relevant ele- 
ments around the HBS in a 1.5Mq, Zq star (typical of the 
lithium-rich RGB giants; see CBOO) which is just passing 
through the RGB bump and for which the extra-mixing 
is still inactive. Since the end of the first dredge-up the 
mass of the core has grown and the convective envelope 
has retreated. At the bump the external wing of the ''Be 
peak crosses the molecular weight discontinuity before the 
""^•^C peak does. When it is connected with the convective 
envelope by the extra-mixing process, ^Be produced via 
^He(a, 7) ""Be starts to diffuse outwards. Depending on the 
mixing timescale r(dif), the CF71 mechanism may or may 
not operate in this radiative region. According to DWOO 
transport coefficients Dr ~ lO^^cm^sec"^ are needed to 
increase the surface lithium abundance. This can be eas- 
ily understood by looking at Fig. ^ where the lifetimes of 
the various elements are compared with T(dif) computed 
for this value of Dr. In particular, we then have T(dif) 
~ T(^Be-|-e~) everywhere below the convective envelope. 
Higher Dr are thus necessary for the freshly synthetized 
^Li to survive, ^Be being transported rapidly to cooler re- 
gions before it decays. In other words, if we assume that 
the extra-mixing is related to rotation and starts with 
the low values obtained by DTOO (10® to lO^Ocm^sec-^), 
then at the beginning of this phase r(dif) > r(''Be-|-e~) : 
''Be diffuses outwards but decays in regions where ''Li is 
rapidly destroyed by proton capture. 

In order to test the reaction of the stellar structure in 
such a situation, we computed an evolutionary sequence 
for a I.5M0, Zq model assuming that the mixing starts at 
the bump and using a parametrized Dr that we take first 
of the order of lO^cm^sec"^ between the base of the con- 
vective envelope and the external wing of the ''Be peak. 
Details of the stellar code and of the input physics will be 
given in Palacios et al. (2001). The energy generation rates 
for the main nuclear reactions are shown in Fig. Before 
the mixing starts ^'*N(p,7)^^0 dominates inside the HBS. 




Fig. 1. 1.5Mq, Zq star at the RGB bump - (top) Lifetime 
of ''Li, 'Be, ^"B and "^He, and convective and diffusive time 
scales (bold lines; T(dif) is shown for a transport coefficient 
constant and equal to lO^^cm^sec^^). The transition be- 
tween the convective enveloppe and the radiative zone is 
indicated by the arrows, (bottom) Abundance profiles of 
the relevant elements (in mass fraction) in the same region 

Then as evolution proceeds with the initially low Dr, ''Be 
diffuses outwards and decays in the radiative layer where 
'Li is rapidly destroyed by proton capture. In this thin 
region the energy released by 'Li(p,a)a increases signifi- 
cantly. This reaction becomes actually the dominant one. 
In addition to the HBS the star has now a very thin lithium 
burning shell (LiBS). In the region where the mixing acts 
the contributions of 'Be(e~, v^lA and ^He('*He,7)''Be also 
increase though to a lesser extent. This change in the en- 
ergy production leads to an increase of the stellar lumi- 
nosity and radius which are rapidly shifted from ~57 to 
~78Lo and from '--^13 to ^16. SR© respectively in the case 
presented here. Consequently the mass loss rate (given by 
Reimers' 1975 approximation in our computations) is al- 
most doubled. One thus expects to see the formation of a 
dust shell around the star as suggested by the IRAS fluxes 
around some of the lithium-rich objects (de la Reza et al. 
1996, 1997). 

In Zahn and Maeder developments for the transport 
of matter and angular momentum, the effective diffusivity 
is proportional to [r U(r)]^, the vertical component of the 
meridional velocity being itself proportional to [Ef^-I-E^]. 
En and E^ depend respectively on the velocity profile and 
on the chemical inhomogeneities along isobars (see e.g. 
Maeder & Zahn 1998 for their complete expressions). But 
what is important for the present discussion is their de- 
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Fig. 2. Energy generation rate [crg.g~-'^.sec~'^] for various 
nuclear reactions in the l.SM©, Z0 star at the RGB bump. 
The H and "^He abundance profiles (in mass fraction) arc 
also shown - (top) Model at the bump, just before the 
start of the extra- mixing, (bottom) Model with mixing. 
enuc(^Li+p) is now the main contributor to the energy re- 
leased in the external wing of the hydrogen burning shell. 
enuc(3He-|-''He) is also significantly increased as well as 
enuc(^Be-|-e~) in the more external layers 

pendence on the local rates of nuclear and gravothermal 

energies. Here we thus propose that the apparition of the 
LiBS and the significant release of enucC Li+p) '^^ highly 
enhance the meridional velocity and the associated trans- 
port coefficient as required for the surface lithium abun- 
dance to increase. 

3. Thermal instabilities generated by lithium 
burning 

In order to mimic this situation wc raised the value of 

Dr (up to 1012-5) 

as 6nuc('''l-'i~t~p) increased. The mixing 
depth was kept constant. We are aware of the crudeness 
of our approach, and of the necessity to follow in a self- 
consistent way the evolution of the transport of matter 
and of angular momentum since the early phases of evolu- 
tion to get the mixing characteristics (see e.g. Charbonncl 
& Talon 1999). However the model used here already gives 
the main trends of the structural response of the star to 
such a mixing. 

In particular, a thermal instability fed by the lithium 
burning, or lithium flash, could develop with the charac- 
teristics shown in Fig. 3. In the thin region where the mix- 
ing acts, enuc(^Li+p) strongly increases (up to ~ 3.2 x 10® 
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Fig. 3. Characteristics of the lithium flash in the 1.5Mq, 
Z0 model in terms of time variations since the beginning 
of the bump (to) of the density (upper panel) , temperature 
(middle panel) and nuclear luminosity (lower panel) at the 
base of the LiBS 

erg g^i sec^i in the model presented here while in a stan- 
dard HBS enuc(CN) is of the order of 2 x 10^ erg g-^ 
sec"'^ at this evolutionary stage). This leads to a fast in- 
crease of the local luminosity. While the LiBS tries to get 
rid of the nuclear energy excess its density drops. How- 
ever its expansion is not sufficient to lift the overlaying 
regions and the corresponding pressure decrease is rather 
small. The temperature thus keeps increasing, as well as 
enuc(''Li+p). Consequently in the thin LiBS the heating 
timescale (cpT/enuc(''Li-|-p)) becomes shorter than the 
timescale of heat diffusion (cpT M^Kp^/{l - /3)cP). The 
thermal instability is then strong enough to require a con- 
vective transport of the energy. As a consequence ^Be can 
be instantaneously transported in the external convective 
envelope. This results in a sudden increase of the surface 
lithium abundance. In the case presented here the lithium 
mass fraction in the convective envelope was increased by 
a factor 3 compared to its post dredge-up value. This 
result is of course dependent on the apparition and on 
the strength of the lithium flash, which itself depends on 
the way the transport coefficient evolves as enuc(''Li+p) 
is dumped in the thin layer between the Be peak and 
the convective envelope. Quantitative conclusions require 
a self-consistent treatment of the whole process. 

A response of Dr to enuc(^Li-l-p) different from the 
one simulated here could lead to a situation where both 
the stellar luminosity sustained by enuc(^Li+p) and the 
surface lithium abundance increase more drastically be- 
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fore the lithium flash. In this ease lithium-rich stars with 
very high luminosity eould be produeed during a very 
short phase. This could explain the rare occurrence of very 
bright lithium-rich giants such as the two stars discovered 
in the globular clusters M3 (Kraft et al. 1999) and NGC 
362 (Smith et al. 1999) which stand at much higher lumi- 
nosity. 

At this stage one can however speculate on the way 
the lithium-rich phase ends. Indeed when the convective 
instability develops, it erases the mean molecular weight 
gradient which was still protecting the deeper layers from 
the mixing. The transport process is then free to proceed 
closer to the core. Once it reaches the regions where -"^^C 
is nuclear ly converted into ^^C (as is observed to occur in 
low-mass stars after the bump) the surface material is ex- 
posed to temperatures higher than lithium can withstand. 
The freshly synthesized lithium is then destroyed as the 
surface carbon isotopic ratio decreases and the lithium- 
rich phase ends. This is consistent with the available data 
on ^^C/^^C in the lithium-rich RGB stars : those with 
the highest lithium abundance still have a standard (i.e. 
post-dredge-up) carbon isotopic value, but no star retains 
its peak lithium abundance once its carbon isotopic value 
dips (CBOO). 

In order to estimate the timescalc of our process let us 
define the beginning and the end of the lithium-rich phase 
as respectively the moments where the surface lithium 
abundance increases above its post-dredge up value and 
where the carbon isotopic ratio decreases below its post- 
dredge up value. While the duration of the lithium flash 
itself is quite short in the model presented here (~ 2 x 10^ 
yrs; see Fig. 3), the hthium-rich phase so-defined lasts for 
~ 60% (i.e. ^ 6 X 10^ yrs) of the classical bump duration. 
This has to be compared with the number of stars located 
in the bump region and which still present a post-dredge 
up carbon isotopic ratio. Among the sample of stars with 
Hipparcos parallaxes and carbon isotopic ratio determina- 
tions used in Charbonnel & do Nascimento (1998), ~ 35% 
of the 49 objects located close to the bump do satisfy this 
criterion. We consider that this comparison is very encour- 
aging in view of the simplifications we made in the present 
model and of the limited observational sample. This pre- 
liminary estimation indicates that many lithium-rich stars 
should be discovered by observational programs focussing 
on the bump region. 

4. Conclusions 

We propose a scenario to explain the lithium-rich phase at 
the RGB bump in the framework of rotation-induced mix- 
ing. When the external wing of the ''Be is connected to the 
convective envelope the extra-mixing timescale is at first 
such that the transported ^Be decays in the region where 
^Li is destroyed by p-capture. A lithium burning shell ap- 
pears. The meridional circulation and the corresponding 
transport coefficient then increase due to the enuc(^Li-l-p) 



burst. ''Li(p,Q:)Q! becomes the dominant reaction leading 
to an increase of the local temperature and of both the lo- 
cal and total luminosities. A lithium flash terminates this 
phase. As the convective instability develops it erases the 
molecular weight gradient and the mixing is free to pro- 
ceed deeper. Both the lithium abundance and the carbon 
isotopic ratio then drop at the stellar surface. During the 
whole sequence the temporary increase of the stellar lumi- 
nosity causes an enhanced mass loss rate which naturally 
accounts for the dust shell suggested by the far-infrared 
color excesses measured for some lithium-rich objects from 
IRAS fluxes. However the contribution of these stars to the 
lithium enrichment of the Galaxy should be very modest. 

The treatment of the transport coefficients used in the 
present simulation is rudimentary and serves only to out- 
line our scenario. In particular the apparition and the 
strength of the lithium fiash depend on the way the trans- 
port coefficient evolves as enuc('^Li-l-p) is dumped in the 
thin layer between the Be peak and the convective enve- 
lope. This has to be followed in a model where the struc- 
ture and the transport of angular momentum and of the 
chemical species are treated sclf-consistently. This will be 
presented in a forthcoming paper where consequences on 
the luminosity function will also be investigated (Palacios 
et al. 2001). 
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